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Abstract 

We consider the calculation of n-point multigluon tree amplitudes with a pair of massive fermions 
in QCD. We give the explicit transformation rules of this kind of massive fermion-pair amplitudes 
with respect to different reference momenta and check the correctness of them by SUSY Ward 
identities. Using these rules and onshell BCFW recursion relation, we calculate the analytic results 
of several n-point multigluon amplitudes. 
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I. INTRODUCTION 



Scattering amplitudes are important from both theoretical and experimental points of 
view. Traditionally we use Feynman diagrams to calculate the scattering amplitudes in field 
theories. During the past several years, motivated by string theory l|, new efficient methods 
for tree level amplitudes have been suggested. The offshell CSW rule ^ suggests that the 
tree amplitudes in gauge theories can be constructed by offshell continued MHV vertices j^. 
Then Britto, , Cachazo, Feng and Witten (BCFW) ^-[6| gave an onshell recursion relation 
where the higher point amplitudes can be obtained by lower point onshell amplitudes whose 
momenta are properly shifted to complex space. 



The BCFW recursion re 



15 



massless gauge theories j?- 



ations have been extensively used to calculate tree level in various 



14| . Its extension to gravity amplitudes has also been considered 



23] . Those applications are all related to massless external particles. But in fact, 



all the matter particles and weak bosons in standard model are massive. Furthermore, 
massive amplitudes are more important in higher energy physics experiments such as LHC 
experiments for processes with top quarks, Higgs particles and possible supersymmetric 



partic 
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es. So it is important to consider the amplitudes with massive external particles. In 
27| . the amplitudes with one external massive gauge bosons or Higgs bosons have been 
discussed. Multigluon amplitudes with pairs of massive scalars or quarks have been studied 
in 
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36] . An excellent and compact expression for multigluon amplitudes with a pair of 



massive scalars or quarks and any number of plus helicity gluons has been found in |38| 



by using an off-shell recursive methods 



39] and the BCFW relations. In j40], the authors 



use some supersymmetric Ward identities to relate a compact expression for multigluon 
helicity amplitudes involving a pair of massive quarks to amplitudes with massive scalars. 
A thorough discussion of BCFW onshell recursion relation for amplitudes with massive 
external particles are given in |32] where the shifted momenta can be massless or massive. 
The authors also use BCFW recursion relation to obtain compact expressions for multigluon 
amplitudes involving a pair of massive quarks or scalars and one minus gluon helicity adjacent 
to fermions. Multigluon amplitudes with a pair of massive scalar and one minus helicity 
gluon not adjacent to scalars are calculated by a different way 



with a pair of massive quarks are calculated in 



30] . All helicity amplitudes 



37] up to six external particles. 



In this paper, we use the notation and convention in [32!| and explore the calculation of 
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several n-point helicity amplitudes with a pair of massive quarks and one minus helicity gluon 
(fermion-pair amplitudes). The massive amplitudes depend on reference momenta which 
define the helicity of the massive fermions. One can relate the fermion-pair amplitudes with 
respect to different reference momenta. We first give the explicit form of the transformation 
rules and check the correctness of them by SUSY Ward identities. Then by shifting the 
momenta of a massive particle and a gluon and using the BCFW recursion relation, we get 
the analytic expressions of several ra-point fermion-pair amplitudes. In section 2, we review 
the spinor formalism of massive fermions and obtain the transformation rules for fermion- 
pair amplitudes defined on different reference momenta. In section 3, we calculate several 
multigluon amplitudes with a pair of fermions. Section 4 is for the summary. 

II. SPINOR FORMALISM OF MASSIVE FERMIONS 

We briefly give our notation and convention of spinor helicity formalism. For massless 
fermions, particles and antiparticles both have definite helicities. Their corresponding spinor 



states are u{p, ±),v{p, ±), which can be denoted as follows |32|: 



u{p, ±) = \pt) , v{p, ±) = \pt) ■ (1) 
For the conjugate states, similar notations are 

u{p,±) = {±p\,v{p,±) = {±p\. (2) 

A massless momentum q'^ can be written in spinor form 

= li-qiriq-) = l{q\r\q] = l{+q\r\q+) = ^[gbl?)- (3) 

The scalar product of two massless momenta p, q can be obtained as 

2p ■ q = {pq)[qp]. (4) 

For a massive momentum = m^, one can always split it into two massless momenta 
by introducing a reference massless momentum q, 

where {K^ f = 0. 



Massive fermions are not helicity eigenstates. Their helicities are frame dependent and we 
can introduce a null reference momentum q to define their helicity states. In this formalism, 
the massive fermions and anti-fermions states with momentum = rn? are 

""^^'^^ " {v^T\q±) 

u{p,±) = + (8) 

t;(p,±) = (gTl^^^^^^(^-m). (9) 
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In the above, p = p" + -^q- We can also define massless states |p ±) = jp^^^^p\qT) and 
rewrite the massive fermion states as 

«(p,±) = |p^T) + ^^r^l^±)- (10) 

We can obtain similar forms for all other states. It is easy to see that we have a smooth 
massless limit. 

Because the massive fermionic helicity states depend on reference momenta, then the 
amplitudes with massive external fermions should also depend on reference momenta. But 
we can relate fermionic spinor states with one reference momentum to these with another 
reference momentum|32| . Let q, q be two light-like reference momenta, we have the following 
relation between the spinor states corresponding to the two reference momenta. 



(11) 



where 

{qp^)[p^q\'' {qP^){p^q)^ [qP'][p^q]' [qP^]{p^q) 

In the above equations, p^{p^) means splitting p with respect to q{q). Similar relations can 
be obtained for other massive fermionic states. 

In this paper, we consider n-point amplitudes with gluons and a pair of massive fermions 
y4(lQ, 2, 3, ■ ■ ■ , — 1,?^q). In the calculation, we choose the same reference momentum 
for the external massive fermions. An amplitude with reference momentum q for the ex- 
ternal fermions is denoted by Aq{lQ, 2, 3, ■ ■ ■ , — 1, nq) . In contrast to amplitudes with 
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massless fermions, there are both hehcity-conserving and hehcity-flipping amphtudes for 
massive amphtudes. With the same external gluons, there are four different helicity am- 



plitudes Aq{lQ, 2, 3, ■ ■ ■ , n 



1 n- 



2, 3, ■ ■ ■ , n — l,n. 



^gi^Q' 2, 3, ■ ■ ■ , n 



Aq{lQ, 2, 3, ■ ■ ■ , n — 1, n^). Similar to the relation flTT]) . we can get relations between ampli- 



tudes with respect to different reference momenta. 



4+" 

A+- 



\ 



( 



\ 



12 



'-'ll'-'ll '-"ll'-"12 '-'12'-"11 '-'12'-': 
"-"ll'-"21 '--'ll"--'22 '-^12'-^21 '-"12'-"22 
"-'2l"-'ll "-"2l"-"12 "-'22"-"ll "-'22"-'12 



\ '-^2l'-^21 "--'2l"--'22 '-^22"-"21 '-"22'-^22 / 



Aq 

V ) 



(13) 



In the above, A^^ is shorthand notation of y4g(lQ, 



and all the amplitudes have the 



same external gluons. As in equation f[T^ . the elements of the matrix are 



(gl^Jg] 



12 



m{qq) 
m{qq) 



' '-^21 



[qpi][p\qy 



\¥n) \piq] ' {¥n) (p^q) ' [qpn] Ip^q] ' [qpn] ipiq) 

These relations between massive amplitudes are important for our calculations in the next 
section and the correctness of it will also be checked there. 



m[qq\ 



[qpi]{p\q) 



(14) 



III. CALCULATION OF FERMION-PAIR AMPLITUDES 



In this section, we use the onshell BCFW recursion relations to calculate several n-point 
amplitudes with a pair of massive fermions and n — 2 gluons. First we list several excellent 
and useful results about amplitudes with a pair of massive particles, which are the building 
blocks of our calculation. 

An excellent and compact expression for the amplitude of a massive complex scalar- 
antiscalar pair and any number of positive helicity gluons is obtained in 38] : 



where 



(n-l)+,n,) = 2-/2-1 



im 



2/1,22/1,3 ■ ■ ■ yi,n-2 (23) (34) ...{n-2,n-iy 



(15) 



Pi,fc 



Pi +P2 H hPfc, 



yi,k = (Pi + P2 H h Pk f - 



(16) 
(17) 



In supersymmetric theories, SUSY transformation can relate a bosonic particle state 
to a fermionic particle state. And there are also relations between amplitudes related by 
SUSY transformations. The SUSY transformations of helicity states have been discussed in 
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42] . In massless SUSY QCD theory, these SUSY transformations have been applied to 



the helicity amplitudes 
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46 1 . I n massive SUSY QCD, similar transformations of helicity 



states have been derived in j40| . Using SUSY transformations, some relations between 
amplitudes with different external helicity particles can also be obtained. These relations 
are the so-called SUSY Ward identities j40|. As already mentioned, in this paper, we choose 
the same reference momentum for the external massive fermions of an amplitude. In this 
case, amplitudes of a massive fermion pair and all plus helicity gluons have the following 
relations with amplitudes where the fermion pair is replaced by massive scalar pair, 

to) 
{p\q) 



A,(lJ,2+,---,(n-l)+,n.) = f^A(l+,2+,---,(n-l)+,n^), 



MIq, 2+, ■ ■ ■ , (n - 1)+, n-^) = ^^(1^, 2+, ■ ■ ■ , - 1)+, n^). (18) 

It is noted that all the SUSY Ward identities as the above are derived by using supersym- 
metry and these relations should hold to any perturbative orders. At the tree level some 
relations can be applied to non-supersymmetric theory. This results from the fact that at 
tree level the quark-gluon amplitudes are the same in both SUSY and non-SUSY theories 
and at loop level there are contributions from SUSY particles to amplitudes. The helicity 
flipping amplitude with all plus external particles vanish, ^1^(1^, 2+, ■ ■ ■ , (n — l)"*", n^) = 0. 
Plugging eq.flTSl) into eq. llTSl) . we can get the basic building blocks for the calculation of 
massive fermion-pair amplitudes. 

There are also similar relations between fermion-pair amplitudes and scalar-pair ampli- 
tudes with one minus helicity gluon 

where Ij = p\\r=p^ is the projection of pi when we choose pj as reference momentum. One 
can note that there is not similar relations for helicity flipping amplitudes when there are 
one minus helicity gluon. 
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A,(l+, 2+, ■ ■ ■ , ■ ■ ■ , (n - 1) + , n-^) = y^A{ll 2+, ■ ■ • , j-, ■ • • , (n - 1)+, n^), 
A,{1q, 2+, . . . ,r , ■ ■ ■ , (n - 1)+, nj) = -A^A(1^, 2+, ■ ■ ■ , j", ■ ■ ■ , (n - 1)+, n±), (19) 



Before we precede to the calculation, here we use the SUSY Ward identities (fT8|) to check 
the correctness of the transformation rules (fT3|) . We know that Aq{lQ, 2"*", ■ ■ ■ , (n— 1)+, Uq) = 
is correct for any reference momentum q. Then using the transformation rules (fT3|) . we 
can obtain 

= iCl,C-J^^ - Cl2C-J^^ + (20) 

where ^(0) = 2+, ■ ■ ■ , (ra — 1)+, n^) . The coefficient of A{(f)) can be calculated as 

J*?] m{qq) {nq) {qI^Jq] m{qq) {lq)_^ m {qq) m {qq) (In) 



ql) [Iq] {qn) {nq) (Ig) {qn) [nq] /qA (Ig) {nq) /qA (Ig) {qn) {nq) m 



{ql) m {qq) {nq) {qn) m {qq) (Ig) ^ {In) {qq) m 



qlj {qn) {nq) {Iq) {qn) (^qlj (ig) {nq) (^qlj (ig) {qn) {nq) 
m {qq) 



qlj {qn) {nq) {Iq) 

= 0. (21) 

So from the transformation rules, we can obtain the correct results Ag{(j)) =0. In a similar 
way, one can check the correctness of the transformation rules (fT3|) from other SUSY Ward 
identities. 

From the SUSY Ward identities, we have obtained the fermion-pair amplitudes with all 
plus helicity gluons. Then we can use onshell BCFW recursion relation to get fermion- 
pair amplitudes with other gluon helicity configurations. In 32|, it has been proved that 



by choosing a proper momenta shift, we can always use onshell BCFW recursion relation 
to calculate fermion-pair amplitudes. In the following, we will calculate several concrete 
multigluon amphtudes. 

The amplitude A{1q, 2^*", ■ ■ ■ , (n— 1)~, n^) can be calculated by shifting momenta 
In spinor formalism, it is 

\n) — > \n) + z\n — 1), 
\n-l] 1^ - 1] - z\n], (22) 



where 



l^±) = |Pn±), (23) 



and 



Pn=Pn + 



m 



-Pn-l- 



2p„-l ■ Pn 

Then using BCFW recursion relation, we can get a compact result as follows, 



(24) 



A„_i(l+,2+,---,(n-l)-,n^) 



n-2 



k=2 

^2n/2-l 



^'^1 ,-Pk,n-l) 



{n, n — 1) 



(1, n - 1) (23) . . . (n - 2, n - 1) S Al^ " 1) 

m2 (fc - 1, k) [2| nS(2/i,, - 1)' 



(25) 



we 



where 5k^2 = 1 — 5fc,2, and when k = 3, 11^=3 (■ ■ ") = 1- From the SUSY Ward identities, if 
miss the factor , we get the corresponding multigluon amplitude with massive scalar. 

It has a more compact form than the one obtained in 30| because we use more compact 
amplitudes eq. (fT5|) as building blocks. But we can check some lower point amplitudes with 
others. The four point scalar amplitude is 



A(l+,2+,3-,4^ 
The five point amplitude is 

A(l+,2+,3+,4-,5^) = t2'/' 

+ ^2^/2 
= ^2^/2 



i2 



(23)pi,3(2|^23^4l3) Phvi 



(26) 



(23) (34)p|4(2|^^^^ 14) 



(4|^3,A|4)' [2|^3,4|4) 



(34)pi,4(3|^3_,^j4) t/i,2(2|^3^,^j4) 

(4|^i^2j4)^ , ,.,3/2 rn^mm^^^ (27) 



(23) {34) pl,{2\p^pjA) 



[34]?/i2|/i^ 



28 



5f 3,4 



2) 



30] up to overall 



These results are the same as the ones obtained from other ways 
conventional coefficients. 

Then we calculate another four multigluon amplitudes with a massive fermion-pair. For 
amplitude 

A(l+,2+,3-,---,(n-l)+,nt), 
we shift the momenta pi and p3. 



|1) ^ |1) + ^|3), 
|3] ^ |3]-z|l]. 



(28) 
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The amplitude A{1q, 2"*", 3 , ■ ■ ■ , (n — l)"*", Uq) can be decomposed as 
As{l^, 2+, 3-, ■ ■ ■ , (n - 1)+, nj) = ^3(1$, 2+, Pstfc, " " " , ^J)^^(-^3:fc, 3", ■ ■ ■ , A;+) 

A:=4 ^3,fc 
n-1 ^ 

+ E ^3(1^, 2+, At;, ■ ■ ■ , n^)—A{-P,-,, 2+, 3-, ■ ■ ■ , /+) 

+ A3(i+,2+,-Pf2)^^3(A+2,3-,---,4). (29) 
P12 

It is easy to see the first two terms in the above equation are both zero because 
there are fermion-pair amphtudes with all plus helicity. Let us see the third term 
^3(1^, 2+, -P,-)^^3(A+2, ■ ■ ■ , n^). We already know 

A3(P+,3-,---,nt) = 0. (30) 

Using the transformation of amplitude with respect to different momenta, it is easy to show 

A3(P+,3-,---,^J) = 0. (31) 

Then 

A3(l3,2+,3-,---,(n-l)+,nt)=0. (32) 
Using the same recursive calculation and induction, we can prove 

2+, ■ ■ ■ , ■ ■ ■ , (n - 1)+, nt) = 0. (33) 

This is consistent with the result from SUSY Ward identities 
by the supersymmetry of massive SUSY QCD. 



40|], which are obtained only 



Then we use the same shifting of momenta as eq.( l28|) and calculate the amplitude 
^3(1^, 2+, 3-, ■ ■ ■ , (n - 1)+, n^). We obtain 



^(1^,2+, 3 



[n - ly,n-^ 

2 

m 



(34) 



(13) (34) ■• -(71-2,71-1) 



X 



n-l 

E 



X { 



{k,k + I){3\p^pj3) 



n-2 
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121 

+ (3l^^^3,.(3|^,^^ - 1]) + 4,.-i(3|^3_„_j2]} 

^ .on/2-1 (^3) 1 g (3|^,^.,J3)3 



(13) (23)---(n-2,n-l) ^p|,(3|^^|^2j2)(3|^^^^J0 

, (M + 1)(3|^,, YT-Uvi, - f,f,,_,)\n - 1] ^ 
yi,r- ■yi,n-2{^\f^f^i\l + l) 



+ i2"/2- 



,i (n3) (3||)J2] 
(13) 2/1,2(23) ■■■(«- 2, n-1) 



(3|^i^j3)(3|^^^^|^3j3)^ 

h ((3|^i^2|3K. + (3|^i^3j3)l/i,2)(3|(2/i,2 + ^,_2^3_p|j) 



{hn-. + 5,^n-.m ,^_^....,^^^_,(3|(,,2+^,_,^3_pb + l) 

For multigluon amplitude ^3(1^, 2+, 3~", • • • , (n— 1)"^, n^), we can use the similar recursive 
method to calculate and we obtain 



>l3(lQ,2+,3-,--.,(n-l)+,nt) (35) 



(n3) (34) •• -(71-2,71-1) 



(3|^.^.J3)^ 



.tlPi,.(3|^,|>3j2)(3|^^^3 JA;)((3||i^^3 J3)yi,2 + (3|^^^^^^ 
/fc^„_i(fc,/c + l)(3||)^|)3j3) -2 

^ ^^^^^^;::;(%^ 

_ .^2-1 (13) 1 g (3|^AJ3)^ 



+ 



(7^3) (23) . . . (n - 2, n - 1) PlMfA,i\'^){^\fA,i\^) 

(M + i)(3||^,,nUi(^i.-^A.-i)l"-^] . 

X {5,._i + Wi- ,,,...,,._2(3|^^^,J^ + 1) > 

_ ■W2-i(13) (3|^J2] 

(7l3) 2/1,2(23) •• -(71-2,71-1) 



X 



(3|^i^2l3)(3|^iA,|3) 



2 



n-l 

S ((3|^i^2|3)pL- + (3|^i^3.|3)yi,2)(3|(yi,2 + ^,^,^3^p|j) 



(i,i + i)(3||^3,.nw(yi.^-i^A.JI^-i] , 

+ ^^^^ _ . .yi,„_2(3|(yi,2 +|>,,|>3_pii + 1) 
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Compairing this result with eq. (lM|) . we can see that they are different from each other just 
by a constant coefficient. And this is consistent with the SUSY Ward identities in eq. (fT9l) . 
Just as eq. (p5|) . eq. (lM|) and eq. (l35|l have more compact forms than the ones obtained from 
other ways. 

Then let us check the massless limit of eq. (1341) . Taking m = 0, eq. (15^ becomes 

A3(l+,2+,3-,---,(n-l) + ,n^) (36) 

^ .on/2-1 (^3) 1 (3|^„_An-l|3)' 

(13) (23) . . . (n - 2,n - 1) J>l,n-i(3||^„_,^,^_j2)(3|^„_^|^^^„_Jn - 1) 
, .W2-i(^3) (3|^J2] 

(13) 2/1,2(23) •• -(71-2,72-1) 

(3|^il>j3)(3||>,An-il3)' 



X 



((3||*A|3)pi,„_i + (3|^i^3^„_j3)yi,2)(3|(yi,2 +|*^^2l*3,„-i)l^ " ^) 



_ (^3)3(31) 
(12) ■■■(nl) 

This is exactly the MHV amplitude with a fermion-antifermion pair in massless QCD theory. 

Finally, we calculate the amplitude ^3(1^, 2"^, 3~ ,■ ■ ■ ,{n — 1)+, Uq) . There is no SUSY 
Ward identity to relate this amplitude with the corresponding scalar one. So this kind of 
amplitudes should be calculated directly from lower point amplitudes with a pair of fermions 
using onshell recursive method. The result is 



^3(1q; 2 ,3 , ■ ■ ■ , 



(n- l)+,n^) 



(37) 



^2n/2-l 



m 



(34)---(r2-2,n-l) 

n-l 

X $:((in) 



(3n) 2 

P3,k 



fc=4 



(3|^3,J1] '''>l.(3|^i^3J2)(3|^i^3,J^)((3|^,|>3j3)yi,^ ^ 



n-2 



(5fc^„_i(A;,/c + l)(3|^^^ |3) 



(3\'0 121 



j2"/2-i 



m(23)-- -(71-2,71- 1) 



1=3 



(377) 



(3|^,,|ir''^pi,(3|^,^,j2)(3|^^^^J/) 



X {5i,n-l + 5/^n-l"^^ 



(/,/ + i)(3|^2^^np'+i(i/i,, - 1] 



2/l,/---?/l,n-2(3|^Aj/ + l) 



■} 
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+ 



(34) . . . (n - 2,n - 1) § mp,^,\S)ph + {S\p^pj3)yi,2){S\(yi,2+P,^,P,.)\j) 



where C(j) is 



+ 



m[23](n|^3|2] m(13)^(n|^3|l] m(13)(3|^j3] 

t/i,2[31](3|^^^j3] |/i,3(23)2(3|^^^j3] (23)(3||ij3](2|(/^ + ^3)!^] 

{13){3\f^\2] yi,3{n\pjl]-ni'{n\p^\l] l/i,2(n|^Jl] 



m(23)i/i,2 ' yi,3(3|^^Jl] |/i,2(3|^Jl] + 2p3-Pn(3|^Jl]- 



(13)(n3)(3|^j2] 
^ m(23)|/i,2(3| (1/1,2 +^,^,|^3_P 13)^^=^'^-^^^'^ (3|^J1] 

In principle, for an n-point massive fermion-pair amplitude with definite minus helicity 
gluon j , we can use the above mentioned shifting to get the analytic expression for it from 
amplitudes with one minus helicity gluon nearer to the fermions. In fact, it is difficult for 
doing it by hand except for some special gluon helicity configurations. But because the 
shift of momenta, recursion decomposition and transformation of amplitudes with respect 
to different reference momenta are all systematic procedures, so it is suitable to develop a 
program to do the work. 



IV. SUMMARY 



In this paper, using onshell BCFW recursion relation and shifting the momenta of a mas- 
sive fermion and a gluon, we calculate several tree level n-point amplitudes with a massive 
fermion-antifermion pair and one minus helicity gluon. Amplitudes with massive fermions 
depend on reference momentum for defining the helicity states of massive external fermions, 
so it is more difficult for calculating them. We give the explicit transformation rules for 
amplitudes with respect to different reference momenta, which are important in the analytic 
calculation of massive fermion-pair amplitudes. The correctness of these rules are checked 
by SUSY Ward identities. We use the most compact and excellent results for some special 
gluon helicity configurations, such as all plus helicity, and pure gluon amplitudes as build- 
ing blocks to obtain four n-point massive amplitudes with more complicated gluon helicity 
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configurations. Generally, calculating the analytic results of n-point massive amplitudes by 
hand is difficult. A program need to be developed to calculate the amplitudes. But it is 
still interesting and calculable to use the recursion method to explore the analytic results 
of amplitudes with finite external particles and more fermions, which are important in high 
energy physics experiments. 
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